Noise-induced vortex reversal of self-propelled particles 
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Wc report an interesting phenomenon of noise-induced vortex reversal in a two-dimensional system 
of self-propelled particles (SPP) with soft-core interactions. With the aid of forward flux sampling, 
we analyze the configurations along the reversal pathway and thus identify the mechanism of vortex 
reversal. We find that statistically the reversal exhibits a hierarchical process: those particles at 
the periphery first change their motion directions, and then more inner layers of particles reverse 
later on. Furthermore, we calculate the dependence of the average reversal rate on noise intensity D 
and the number N of SPP. We find that the rate decreases exponentially with the reciprocal of D. 
Interestingly, the rate varies nonmonotonically with TV and a minimal rate exists for an intermediate 
value of N. 

PACS numbers: 05.65.+b, 87.18. Tt, 64.60.-i 
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In recent years, the collective dynamics of self- 
propelled particles (SPP) has been a subject of intense 
research due to the potential implications in biology, 
physics, and engineering (see [1[ for a recent review). 
Examples of SPP are abundant including traffic flow |2j , 
birds flocks [!, H| , insects swarms [El Bl, b acteria colonies 
@, and active granular media |9Ul0j . just to list a 
few. Especially, inspired by the seminal work of Vicsek 
et al. [11], many theoretical and experimental studies re- 
ported that various systems of SPP can exhibit a wealth 
of emergent noncquilibrium patterns like swarming, clus- 
tering and vortex [Ia - [l8| . 



A fascinating phenomenon about SPP is that they can 
abruptly change their collective motion pattern, which 
may be induced by either intrinsic stochasticity such as 
error in communication among SPP or a response to an 
external influence such as a predator. For example, a re- 
cent experiment showed that marching locusts can sud- 
denly switch their direction without any change in the 
external environment [f|. Later, the experimental re- 
sults are further explained theoretically by a mathemati- 
cal modeling, highlighting the nontrivial role of random- 
ness or noise on this transition (l9| . Also, noise- induced 
transitions between translational motion and rotational 
motion of SPP have been observed [2CH23| . However, the 
research on this topic is still in its infancy and deserves 
more investigations. In particular, investigation of the 
mechanisms about these transitions is lacking at present. 
As we know, identifying the underlying mechanisms is 
key to understanding and controlling the collective mo- 
tion of SPP. 

In this paper, we report an interesting phenomenon of 
noise-induced vortex reversal between two different ro- 
tational directions in a two-dimensional model of SPP 
interacting via Morse potential [l7|. Vortex pattern has 



been commonly observed in nature, such as fish, ants 
24 , [25j . Bacillus subtilis [26| . and Dictyostelium cells 
27j . By virtue of a recently developed simulation method 
of rare event, forward flux sampling (FFS) [28j . we an- 
alyze the intermediate configurations along the reversal 
path and compute average reversal rate. We find that 
an important statistical property of the vortex reversal, 
that is, the reversal first starts from peripheral particles 
and then gradually to inner particles, so that almost all 
particles change their rational directions. Furthermore, 
we show that the reversal rate decreases exponentially 
with the inverse of noise intensity. Interestingly, the rate 
varies nonmonotonically with the number of particles and 
a minimal rate exists. 

We consider TV identical SPP in two-dimensional space 
with positions £4 , velocities Vi [i = "L - ■ ■ , N) and unit 
mass. The equations of motion read [13] 
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where the first and second terms on right hand side (rhs) 
of Eq|2] represent self-propelled force and friction force, 
respectively, and the third term is pair interaction among 
SPP given by the generalized Morse potential 
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Here, l a and l r represent the attractive and repulsive 
potential ranges, C a and C r represent their respective 
amplitudes. The last term on rhs of Eq[2] is a stochas- 
tic force of intensity D that are independent of particle 
index and satisfy (&,m(*)) = and (&,m = 
SijS mn S (t — t') with i, j £ 1, • • • , N; m,n € x,y. 

The model exhibits diverse dynamic patterns, such as 
clumps, rings and vortex, for different model parameters. 
Here we set a = 1.0, (3 = 0.5, l a = 2.0, l r = 0.5, C a = 0.6 
and C r = 1.0, which corresponds to the case of vortex. 
To distinguish two different rotational directions of vor- 
tex, we define average angular momentum of particles 
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FIG. 1: (Color online) A long-time evolution of L(t) at TV = 
40 and D — 0.37 shows several vortex reversals. The inset 
depicts two typical stable configurations of SPP: vortex with 
CW rotation (L < 0) and CCW rotation (L > 0). 

L(t) as L(t) = J2iLi Li(t), where Li(t) = Xi(t) x 
is angular momentum of particle i at time t. L < 
and L > indicate that SPP rotate clockwise (CW) and 
counterclockwise (CCW), respectively. In the presence 
of a weak noise, vortex pattern is robust to noise. 

Figure [T] shows a long-time evolution of L(t) at TV = 
40 and D = 0.37, obtained by numerically integrating 
Eqs lll2l using a fourth order Adams-Bashforth method 
allowing particles an infinite range of motion. One 
can observe that noise can induce the sudden transitions 
of SPP between two different rational directions. The in- 
set in Fig[T] depicts two typical configurations of SPP: 
vortex with CW and CCW rotation. However, these 
transitions occur rarely and average waiting time be- 
tween transitions is very long. In this situation, conven- 
tional brute- force simulation becomes highly inefficient. 
To overcome this difficulty, we will use FFS method of 
Allen and coworkers [28[ to compute the rate of vortex 
reversal and evaluate statistical properties of the rever- 
sal path, which is the main purpose of the present work. 
FFS method was designed to study rare events both in 
and out of equilibrium. This method first defines an or- 
der parameter to distinguish between the initial state A 
and the final state B, and then uses a series of interfaces 
to force the system from A to B in a ratchet-like manner. 
Here, it is convenient to select L as the order parameter, 
and consider CW vortex as A and CCW vortex as B 
without loss of generality because these two states are 
equivalent in our model. 

By storing SPP configurations at each interface of FFS 
sampling, one can identify the statistical properties of the 
vortex reversal pathway. All results below are obtained 
by averaging 10 independent FFS samplings. In each 
FFS sampling, 1000 configurations in the center of mass 
coordinates are stored at each interface and are analyzed 
to obtain the statistical properties of the configurations. 
In Fig|3J we show the velocity field of SPP at six different 



FIG. 2: (Color online) The velocity field of SPP at six dif- 
ferent FFS interfaces, corresponding to (a) L — —1.95, (b) 
L = -1.6, (c) L = 0, (d) L = 0.25, (e) L = 0.4, and L = 1.0, 
respectively. The bar in (a) indicates unit length. 



interfaces. Figj2ja) and FigJSJf) show the velocity field 
before and after the reversal, indicating that the system 
is in stable vortex with CW and CCW rotation, respec- 
tively. While Figsl^b-e) show the intermediate processes 
of the reversal. One can clearly observe that along the 
pathway of vortex reversal the velocity field gradually 
changes its sign from the periphery to center of the vor- 
tex. That is to say, vortex reversal first starts from pe- 
ripheral particles and then gradually to inner particles, 
and finally almost all particles change original rational 
directions. 

Further information on pathway of vortex reversal is 
provided by the radial distributions of average angu- 
lar momentum L(r) and average angular velocity w(r), 
where r is the distance to the center of mass. In Fig [31 
we plot L(r) and ui(r) for five different interfaces along 
the pathway of vortex reversal. From the variations of 
L{r) and u)(r) with interfaces one can observe that the 
whole process of the vortex reversal. Before the rever- 
sal (L = —1.95) the values of L(r) and oj{r) are always 
negative, irrespectively of r. That is, the system is in sta- 
ble vortex with CW rotation before the reversal. When 
the reversal happens, for example, for L = the val- 
ues of L{r) and w(r) become positive for r > 1.8, while 
for r < 1.8 they are always negative. With increasing 
L this situation further goes till any values of L(r) and 
uj(r) become positive. Finally, the vortex rotates with 
CCW direction. Therefore, this further validate that the 
statistical property of vortex reversal is that the reversal 
process starts from the periphery of the vortex. Also, 
we calculate the radial distributions of density of parti- 
cles and the absolute value of velocity of particles, and 
find that they do not have significant difference with the 
interfaces L. 

It is informative to analyze the statistical properties 
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FIG. 3: (Color online) Radial distributions of average an- 
gular momentum L(r) and average angular velocity uj(r), for 
five different interfaces along the pathway of vortex reversal. 
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FIG. 4: (Color online) The distributions of angular momen- 
tum Li near the critical configurations for these particles with 
r £ [0.6, 1.0] (squares), r £ [1.6, 2.0] (circles), and r £ [2.6, 3.0] 
(triangles) . 



of critical configurations of SPP. Similar to definition in 
previous studies [3(j, the critical nucleus is determined 
by the committor probability -Pb(*) = 0.5, where Pb(i) 
is the probability of reaching B state before returning 
to A state starting from the FFS interface i, which can 
be computed by FFS sampling. This shows that when 
the system initially locates at the critical configurations 
there is equal probability of returning to CW vortex or 
CCW vortex. We find that L ~ 0.08 at the critical con- 
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FIG. 5: (Color online) The natural logarithm of rate of vortex 
reversal In R as a function of the inverse of noise intensity 1/D 
for different SPP number N (a) and as a function of N for 
different D (b). The dashed lines in (a) are plotted by linear 
fitting. 



figurations, that is, the average angular momentum of all 
particles is slightly larger than zero at the critical config- 
urations. By varying noise intensity D and the number 
of particles N, the value is not almost changed. Further- 
more, we analyze the radial distributions of angular mo- 
mentum Li of particles near the critical configurations. 
In Fig. 4, we show that the distributions of Li of par- 
ticles at r G [0.6,1.0] (squares), [1.6,2.0] (circles), and 
[2.6,3.0] (triangles). Interestingly, all these distributions 
are bimodal that are independent of r. The two peaks 
always correspond to Li < and Li > 0, respectively, 
separated by the lowest value of distributions located 
at about Lj ~ 0. But the loci of two peaks and their 
relative heights vary with r. For particles at the inner 
layer of vortex, the absolute values of is relatively 
low, meaning that the rotational property of these parti- 
cles is weaken and thus become more disordered. While 
the other particles are of obvious ration since the values 
of Li are comparable with ones at stable vortex. The 
first peak for Li < is higher at the inner layer of vor- 
tex, while the second peak for L t < is higher at the 
outer layer of vortex. At the center layer of vortex, the 
heights of the two peaks are nearly the same. The results 
give the properties of critical configurations and thereby 
further illustrate the mechanism of vortex reversal. 

Another key question in vortex reversal concerns the 
dependence of average rate R of the reversal on noise 
intensity D and the number of particles N. The FFS 
method can give R as the production of the flux across 
the first interface and the conditional probabilities of 
reaching the last interface without returning the first in- 
terface. In FigGJa) we plot the natural logarithm of the 
rate In R as a function of the inverse of D for N. The 
results show that R decreases exponentially with 1/D. 
It should be pointed out that in our simulation D is nei- 
ther too large nor too small. We find that if D is larger 
than a critical value D c ~ 0.44, vortex pattern will be 



4 



destroyed, and thus vertex reversal makes no sense. If 
D is too small, R is very low so that simulation will be 
time-consuming. In FigJSJb) we show that In R as a func- 
tion of N for different D. Interestingly, we find that R 
varies nonmonotonically with N . There exists a mini- 
mal value of R at N = 48 for low D and at N = 50 
for high D. From FigJSJa), one may speculate that there 
seems to be an effective nonequilibrium potential for de- 
scribing the collective transition in rotational directions. 
The potential barrier between CW vortex and CCW vor- 
tex is fixed if N is unchanged, such that the transition 
rate follow classical Kramers' law. On the other hand, 
the nonmonotonic dependence on R ~ N implies depen- 
dence of the effective nonequilibrium potential on N is 
nontrivial if exists. However, understanding these results 
from the view of theoretical analysis, if not infeasible, is 
at least a complex task at present. 

In summary, using a two-dimensional model of SPP in- 
teracting via a soft-core potential, we have investigated 
the mechanism of noise-induced the changes of vortex 
pattern in rational direction. By virtue of FFS method 
we analyze the statistical property and compute the rate 



of the reversal. We find that the reversal process is hier- 
archical: the process initially inspired by the peripheral 
particles, and those particles gradually drive more inner 
layers of particles into reverse motion directions. On the 
other hand, we show that the rate of the reversal de- 
creases exponentially with the inverse of noise intensity. 
Interestingly, the reversal rate depends nonmonotonically 
on the number of SPP and a minimal rate exists at a mod- 
erate number of particles. Our findings may provide us 
some new understanding on the transitions of collective 
patterns of SPP. 



Acknowledgments 

This work is supported by NSFC (Grant Nos. 
91027012, 20933006). Z.H. acknowledges support by 
China National Funds for Distinguished Young Scientists 
(Grant No. 21125313). H.C. acknowledges support by 
the Doctoral Research Foundation of Anhui University 
(Grant No. KJ2012B042). 



[1] T. Vicsek and A. Zafiris, e-print: larXiv:1010.50T7l f2010h 
[2] D. Helbing, Rev. Mod. Phys. 73, 1067 (2001). 
[3] A. C. et ai, Proc. Natl. Acad. Sci. U.S.A. 107, 11865 
(2010). 

[4] K. Bhattacharya and T. Vicsek, New J. Phys. 12, 093019 
(2010). 

[5] J. B. et al, Science 312, 1402 (2006). 

[6] P. Romanczuk, I. D. Couzin, and L. Schimansky-Geier, 

Phys. Rev. Lett. 102, 010602 (2009). 
[7] L. Angelani, R. Di Leonardo, and G. Ruocco, Phys. Rev. 

Lett. 102, 048104 (2009). 
[8] E.-L. F. H. P. Zhang, Avraham Be'er and H. L. Swinney, 

Proc. Natl. Acad. Sci. U.S.A. 107, 13626 (2010). 
[9] S. R. V. Narayan and N. Menon, Science 317, 105 (2007). 
[10] J. Deseigne, O. Dauchot, and H. Chate, Phys. Rev. Lett. 

105, 098001 (2010). 
[11] T. Vicsek, A. Czirok, E. Ben- Jacob, I. Cohen, and 

O. Shochet, Phys. Rev. Lett. 75, 1226 (1995). 
[12] J. Toner and Y. Tu, Phys. Rev. Lett. 75, 4326 (1995). 
[13] G. Gregoire and H. Chate, Phys. Rev. Lett. 92, 025702 
(2004). 

[14] M. Aldana, V. Dossetti, C. Huepe, V. M. Kenkre, and 
H. Larralde, Phys. Rev. Lett. 98, 095702 (2007). 

[15] H. Chate, F. Ginelli, and R. Montagne, Phys. Rev. Lett. 
96, 180602 (2006). 

[16] H. Levine, W.-J. Rappel, and I. Cohen, Phys. Rev. E 63, 
017101 (2000). 

[17] M. R. D'Orsogna, Y. L. Chuang, A. L. Bertozzi, and L. S. 
Chayes, Phys. Rev. Lett. 96, 104302 (2006). 



[18] F. Peruani, T. Klauss, A. Deutsch, and A. Voss-Boehme, 

Phys. Rev. Lett. 106, 128101 (2011). 
[19] C. A. Yatesa, R. Erbana, C. Escudero, I. D. Couzin, 

J. Buhl, I. G. Kevrekidis, P. K. Mainia, and S. D. J. 

T, Proc. Natl. Acad. Sci. U.S.A. 106, 5464 (2009). 
[20] A. S. Mikhailov and D. H. Zanette, Phys. Rev. E 60, 

4571 (1999). 

[21] U. Erdmann, W. Ebeling, and A. S. Mikhailov, Phys. 

Rev. E 71, 051904 (2005). 
[22] J. M. Allison Kolpas and I. G. Kevrekidis, Proc. Natl. 

Acad. Sci. U.S.A. 104, 5931 (2007). 
[23] J. Strefler, U. Erdmann, and L. Schimansky-Geier, Phys. 

Rev. E 78, 031927 (2008). 
[24] S. V. V. Julia K. Parrish and D. Gmnbaum, Biol. Bull. 

202, 296 (2002). 
[25] D. J. T. Sumpter, Phil. Trans. R. Soc. B 361, 5 (2006). 
[26] A. Czirok, E. Ben- Jacob, I. Cohen, and T. Vicsek, Phys. 

Rev. E 54, 1791 (1996). 
[27] W.-J. Rappel, A. Nicol, A. Sarkissian, H. Levine, and 

W. F. Loomis, Phys. Rev. Lett. 83, 1247 (1999). 
[28] R. J. Allen, P. B. Warren, and P. R. ten Wolde, Phys. 

Rev. Lett. 94, 018104 (2005). 
[29] G. H. Golub and J. M. Ortega, Scientific Computing 

and Differential Equations: An Introduction to Numeri- 
cal Methods (Academic Press, 1992). 
[30] A. C. Pan and D. Chandler, J. Phys. Chem. B 108, 19681 

(2004). 



